Non-perturbative microscopic theory of superconducting fluctuations near a quantum critical point 
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We consider an inhomogeneous anisotropic gap superconductor in the vicinity of the quantum critical point, 
where the transition temperature is suppressed to zero by disorder. Starting with the BCS Hamiltonian, we 
derive the Ginzburg-Landau action for the superconducting order parameter. It is shown that the critical theory 
corresponds to the marginal case in two dimensions and is formally equivalent to the theory of an antiferro- 
magnetic quantum critical point, which is a quantum critical theory with the dynamic critical exponent, z = 2. 
This allows us to use a parquet method to calculate the non-perturbative effect of quantum superconducting 
fluctuations on thermodynamic properties. We derive a general expression for the fluctuation magnetic suscep- 
tibility, which exhibits a crossover from the logarithmic dependence, x(T, H, n) oc In [6n(T, H)], valid beyond 
the Ginzburg region to x(T,H,n) oc ln I/5 [Sn(T, H)] valid in the immediate vicinity of the transition (where Sn 
is the deviation from the critical disorder concentration). We suggest that the obtained non-perturbative results 
describe the low-temperature critical behavior of a variety of diverse superconducting systems, which include 
overdoped high-temperature cuprates, disordered p-wave superconductors, and conventional superconducting 
films with magnetic impurities. 
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The problem of the critical behavior of itinerant electronic 
systems in the vicinity of quantum phase transitions has 
been the subject of intense theoretical investigations in recent 
years III S |3l |4j, |5[] . Most theoretical studies of the problem 
are based on effective field theories, which describe critical 
fluctuations near the transition. The bare Ginzburg-Landau 
coefficients of these models are usually treated as phenomeno- 
logical parameters, which often remain undefined due to the 
lack of a reliable microscopic theory of the transition. Indeed, 
while Fermi liquid theory predicts and classifies possible in- 
trinsic instabilities, there is usually no controlled approach to 
access the transition point on the basis of a fermionic Hamil- 
tonian. The notable exception is a superconducting instability, 
which is well-described by the perturbative Bardeen-Cooper- 
Schrieffer (BCS) theory. A strong magnetic field or disorder 
effects may suppress the superconducting transition tempera- 
ture to zero and therefore lead to a superconducting quantum 
critical point (QCP), which can be studied on the basis of the 
BCS model. 

The critical behavior of a superconducting system near 
the transition is governed by superconducting fluctuation ef- 
fects |6J], which physically are due to uncondensed fluctuating 
Cooper pairs, which co-exist with electronic excitations. A 
perturbative theory of classical superconducting fluctuations 
was developed by Aslamazov and Larkin in 1968 |7]. More 
recently, Larkin and the author [8] considered quantum super- 
conducting fluctuations near the magnetic-field-tuned QCP. In 
both cases, it was found that Gaussian fluctuations strongly 
affect thermodynamics and transport near the critical point. 
Even though the Alsamazov-Larkin theory has been very suc- 
cessful in describing a variety of experiments, it is strictly 
speaking not a critical theory, but a Gaussian perturbation the- 
ory, which assumes that the fluctuating Cooper pairs do not in- 
teract, and applies only far enough from the transition (beyond 
the Ginzburg region). To the best of the author's knowledge, 
there is no physically relevant example of a non-perturbative 



treatment of superconducting fluctuations. 

In this paper, we develop such a non-petrubative micro- 
scopic theory for a disordered anisotropic gap superconduc- 
tor near the disorder-tuned QCP. We derive the corresponding 
critical theory, which is shown to be identical to the Hertz- 
Moriya-Millis theory J2I 0] of an antiferromagnetic QCP in 
two dimensions. We find that the dimensionless bare quar- 
tic coupling (which characterizes the interaction between su- 
perconducting fluctuations) is a small parameter of the or- 
der of the inverse conductance and this quartic term becomes 
marginally irrelevant at the transition. This allows us to per- 
form non-perturbative parquet resummation of the leading 
logarithms and find the exact critical behavior of the magnetic 
susceptibility near the transition. The latter crosses over from 
X(T, n) oc In [n - n c (T)/n] (which is the quantum analogue of 
the Aslamazov-Larkin result valid far enough from the tran- 
sition) to the critical behavior x(T, n) oc ln 1 ^ 5 [n — n c (T)/n], 
which holds within the quantum Ginzburg region. 

Let us consider a disordered electron system with a weak 
attraction in the /-wave channel, described by the Hamiltonian 

ik = Y^tUptp + \ X v (P.P'WWp'#-p'+q + #dis,(D 
p pp'q 

where ^ p = £(p) - /u is the spectrum, 'Kdis represents a dis- 
order potential (which we assume to be due to Poisson dis- 
tributed short-range impurities), the interaction is V(p, p') = 
-Aiui((p)ui(<p), A/ is the interaction constant, and Ui(<f>) de- 
fines the symmetry of the gap. In what follows we assume 
an unusual pairing symmetry (e.g., af-wave), so that TTi = 

jj " u,((f>)d(f>/(2n) = and u] = 1. We suppress spin indices 
throughout the paper. Next, we introduce an anisotropic order 
parameter and allow for its spatial and temporal fluctuations 

A k (r, r) = T £ y ( k > k ')^ ( k ' - \ • k + § ; w «) e ' q '™> 

(2) 
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where t is the Matsubara time and F(r, r;r',r') = 
- (T T t[r(r, r)iWr, t)^ is the Gor'kov's Green's function. Fol- 
lowing Ref. [9], we assume that the symmetry of the gap is 
preserved, Ak(r, r) = A(r, T)M/(</>it), and integrate out the one- 
particle degrees of freedom to obtain the following effective 
action for the order parameter near the transition 
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where v is the density of states at the Fermi line and we 
use a three- vector Q = (a>„,q) and the notation fgf(Q) = 

TYj J d 2 q/(2n) 2 f(ci> n ,q) for brevity (oj„ is the bosonic Mat- 

subara frequency). The general expression for the quartic co- 
efficient B (at an arbitrary temperature and disorder strength) 
was derived by the author in Ref. 19J] . Near the QCR Tt c q — > 
0, it reduces to B = vT 2 (rf - 1/3) [c.f., B = 7^(3)v/(8^ 2 r 2 ) 
in the classical limit, Tt c q » 1], where the overline implies an 
averaging over the directions on the Fermi line. E.g., in a d- 
wave superconductor, we get = 3/2. In Eq. ©, the critical 
scattering time (or disorder concentration, n c ~ t^ 1 ) deter- 
mines the superconducting transition point and is a function 
of temperature and magnetic field. To find the dependence 
on the latter, we can just replace the Cooper-pair momen- 
tum q with the operator ft - [-/V - — A(r)J and evaluate the 
matrix element corresponding to the lowest Landau level in 
Eq. © by replacing Dq 2 with (0|Dtt 2 |0) = leDH/c. The full 
three-dimensional temperature-disorder-magnetic field phase 
diagram is determined by the equation 
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where tfr(z) is the logarithmic derivative of the gamma- 
function. The zero-temperature solution of this equation rep- 
resents a line of quantum critical points. In what follows, we 
consider the low-temperature limit and small magnetic fields. 
In this case, the asymptotic behavior of the critical scattering 
time is 



t c (T, H) - T c0 
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where g = Epr/n » 1 is the dimensionless conductance, co c = 
eH/(mc) is the cyclotron frequency, 1/t c o = nT c o/y is the 
critical scattering rate corresponding to T = and H = 0, 
y ~ 1.781 is the Euler's constant, and T c q oc exp[-l/(/l;v)] is 
the BCS transition temperature in a clean system. 

Now, we consider the ultra-low-temperature regime in 
which quantum rather than thermal fluctuations determine the 
critical behavior of the system, which corresponds to the limit 
HJ,H,t) = [t c (T,H) - t] /t » 2nTr. In this case, we can 
replace the Matsubara sums in Eq. © with the integrals over 
the frequency. It is also convenient to introduce new variables 



r = [t -t c (T,H)] It and p = B/^DtV) and rescale the 
parameters in Eq. (O as follows: kg = cot, k = VOrq, and 
<p = -\/2v/(Dt 2 )A. Using these dimensionless variables, we 
arrive to the following quantum critical action 

SW= ^ J k [r+\k \+k 2 ]\<p(k)\ 2 (6) 

p u r f(h)<p(k 2 )f(h)<p(k x +h-k 2 ), 



where the parameter po = (3u'j - l)/g <K 1 is essentially the 
inverse conductance and thus is small. We emphasize that the 
above action is not a phenomenological effective theory, but a 
microscopic result, which follows from the basic BCS Hamil- 
tonian ([T). The latter contains just one phenomenological 
parameter - the Cooper channel interaction constant, which 
is encoded in a non-universal "clean" transition temperature, 
r c o. The exact mechanism of unconventional superconductiv- 
ity and the exact value of T c o are not important for the prob- 
lem at hand, which deals with disorder-induced suppression 
of superconductivity that happens at the diffusive rather than 
ballistic length-scales. Another important observation is that 
action © is identical to the Hertz-Moriya-Millis theory QUI 
of the two-dimensional antiferromagnetic QCR The latter is 
known to be a marginal theory II 1011 . Since the dimension- 
less coupling constant in Eq. (O is small, the renormalization 
group (RG) treatment of the theory is an asymptotically exact 
approach. 

The upper cut-off in the RG approach is determined by the 
applicability of the diffusion approximation, which is reason- 
able as long as |w max | ~ Dq 2 ^ < 1/r. In terms of the dimen- 
sionless three-momentum k = (ko, k), this condition implies 
l^maxl = A ~ 1. We use this value of the cut-off and the 
usual RG scheme [01, to obtain the scattering amplitude of 
fluctuations p(k\,k2,k^,kji). In principle, the latter depends 
on all external momenta, but within the logarithmic accu- 
racy, we can just consider it to be a function of a single vari- 
able p{k- max {&,}), for which we get the following equa- 
tion Ji M p(k) = po - I J k p 2 (k') (r + |jf | + k' 2 )~ 2 , which 
leads to the "zero-charge" behavior of p(k) 
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Since the bare scattering amplitude is determined by the in- 
verse conductance, it is alluring to interpret this "zero charge" 
result as a flow of the resistivity to zero in the superconduct- 
ing phase. In fact, a naive calculation of the Aslamazov- 
Larkin conductivity diagram J7|] in which the current vertices 
are taken to be independent of frequencies does lead to the 
logarithmic behavior of the corresponding correction to the 
conductivity 111 ill , which would be consistent with the above- 
mentioned interpretation. However, a more careful analysis 
of the Alsamazov-Larkin diagram shows that this standard ap- 
proximation [6] (valid near a classical transition) breaks down 
near a QCP (see also Ref. [8], where this complication was 
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first pointed out) and a full frequency and momentum depen- 
dence of the current vertices is needed to get a correct result. 
To calculate transport properties turns out to be very difficult 
due to the problem of analytical continuation of the propa- 
gators and vertices (which here are very complicated non- 
analytic functions of two complex variables). However, to 
obtain thermodynamic properties is rather straightforward and 
can be done non-perturbatively on the basis of the action ©. 

Below, we consider the fluctuation magnetic susceptibility 
near the quantum phase transition. 
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where the r(//)-dependece is given by Eq. (0, which yields 
(dr/dH) = eT 2 (vl/c 2 ). Now, we follow Larkin and Khmel- 
nitskii [12] and notice that the second derivative of the free 
energy in ([8]l is the exact polarization operator given by 
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where T(k) is the vertex function, which is determined by 
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Eqs. (O, (O, and ( TTOb lead to the magnetic susceptibility per 
unit volume 
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where d is the thickness of the film (or interlayer distance in 
the case of a layered superconductor) and the proximity to 
the transition r(T, H, t) is given by Eq. ©. Eq. ( fTTT i has two 
asymptotic regimes, 



X 



e ( (g/n) |lnr| , if 1 <K | lnr| <k g; 

dhm^ \ 7.358 «^ 5 |lnr| 1/5 , if|lnr|»g, 



(12) 



where the second numerical coefficient corresponds to the d- 
wave case. The former (finite r) behavior is the regime of 
regular Gaussian perturbation theory. The later regime of 
r — > is clearly non-perturbative and becomes asymptot- 
ically exact in the very vicinity of the QCP. Note that the 
fluctuation diamagnetism exceeds the Landau susceptibility 
AfLandau = —e 2 1 {\2ndhmc 2 ) by orders of magnitude even far 
from the transition. We reiterate that there is no contradiction 
here |6[] since Eq. (fTTT ) should be viewed as a correction to the 
perfect diamagnetic response of a superconductor (not to the 
Pauli/Landau terms in a normal metal). We also note here that 
the effect of quantum fluctuations on other thermodynamic 
properties, such as the specific heat, are unremarkable since 
the critical density of disorder depends on the temperature as 
r(T) oc T 2 . 



Now we discuss the crossover between the quantum 
fluctuation and classical fluctuation regimes. Clearly at 
large temperatures, the non-perturbative RG treatment breaks 
down, because the integral over frequency (which makes 
the quantum problem effectively 4-dimensional) is replaced 
with the Matsubara sum and only the co n = term 
plays a role near the transition; thus, we recover the two- 
dimensional model. The crossover between the two behav- 
iors (here we are talking about the linear-response H — > 
magnetic susceptibility) happens around r ~ 2jiTt = 
2y(T/T co ). The general expression for the leading order 
correction to the fluctuation susceptibility is determined by 

X W = - (dr/dH) 2 T £ w „ J d 2 q/(2n) 2 [r + \cj n \r + Dq 2 r]~ 2 , 



which leads to the result 



X W (T) = - 
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The low-temperature quantum limit in Eq. ( fT3l ) reproduces 
the logarithmic asymptotic of Eq. ( fT2l . while the "high- 
temperature" limit r <sc (T/T c q) <k 1 leads to the classical 
Aslamazov-Larkin-type power law divergence 
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where we introduced a scattering time-dependent transition 
temperature [which is the inverse function of t c (T) used ear- 
lier, see Eq. ©]. We emphasize that Eq. ( TT4T > corresponds to 
the limit (T/T c q) <k 1. However, at higher temperatures and 
in particular near T c o (clean limit), the familiar Aslamazov- 
Larkin power law and all parameters are preserved and only 
the numerical coefficient changes. 

Fig. 1 summarizes the behavior of fluctuations in the vicin- 
ity of the disorder-induced superconductor-metal transition 
in an anisotropic gap superconductor. The solid black line 
represents the phase boundary between the superconducting 
and metallic phases. The hatched area is the Ginzburg re- 
gion, where the fluctuations interact strongly and Aslamazov- 
Larkin theory breaks down. Interestingly, the width of the 
quantum Ginzburg region, G1q = exp(-g), is much smaller 
than that of the classical Ginzburg region, Gic = l/g- The 
two dashed lines in Fig. 1 separate the classical and quan- 
tum fluctuation regimes. The classical regime is effectively 
a two-dimensional theory (d^ = 2) in which the leading or- 
der perturbative correction diverges as a power law, Xr oc 
-(T — T c y l . The effective dimensionality in the quantum 
regime is d 1 ^ = 4 and the leading order correction to sus- 



ceptibility (and likely conductivity 111 111 ) is logarithmic, Xq* oc 
In (n - « c ). The region in between the two dashed lines rep- 
resents a crossover between the quantum and classical behav- 
ior JHIi, 13] and in some sense describes a crossover between 



the effective dimensionality — 4 in the former to the ef- 
fective dimensionality cF s = 2 in the latter. The leading order 
correction to thermodynamics in this regime is described by 
the non-linear function in Eq. ( fT3l , which we found to be a 
smooth function with no remarkable properties. The exact 
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FIG. 1: (Color online) This figure illustrates the phase diagram (tran- 
sition temperature vs. concentration of disorder) and various fluctu- 
ation regimes. The hatched region is the regime of strong supercon- 
ducting fluctuations. The dashed lines determine the quantum, clas- 
sical, and quantum-to-classical crossover regions in the fluctuation 
thermodynamics. See text for a more detailed discussion. 

non-perturbative behavior in the classical and crossover re- 
gions is unknown, but some insight can be obtained via an 
e-expansion J3, ■ In the quantum Ginzburg region, the RG 
approach is asymptotically exact and leads to the susceptibil- 
ity, which behaves as^Q oc ln 1 ^ 5 (n - n c ). 

We note that in a disordered system, the latter non- 
perturbative behavior may be smeared out by the Griffiths ef- 
fects due to mesoscopic fluctuations of disorder. There always 
exist superconducting puddles, which appear in the local ar- 
eas, where the concentration of impurities is smaller than the 
average H5I1 . At low temperatures, these puddles are 

connected via Josephson tunneling, which is long-range in 
real space 7(r (/ ) oc |r,- - r,-|~ 2 . A global transition in this ran- 
dom network is governed by phase fluctuations 11411 . The cor- 
responding theory is however not a pure random .XT-model. 
First, in a af-wave superconductor (or any other supercon- 
ductor where the gap vanishes), there always exist gapless 
quasiparticles with a small density of states, which lead to 
an unusual intrinsic dissipation and non-local in time Joseph- 
son coupling even in the absence of a normal component (the 
long-time asymptote of both the dissipation and Josephson 
kernels behave as lfl6ll a(r) oc J(t) oc |rp 3 ). For supercon- 
ducting puddles surrounded by a normal metal with a large 
electronic density of states, the dissipation kernel acquires the 
usual Caldeira-Leggett form a(r) oc |t|~ 2 , while the Joseph- 
son term remains long-range in time J(t) oc |t|~ 3 . Within a 
mean field approach, one can estimate [11] that the QCP © 
is shifted by 6r ~ g~ l and therefore the quantum Ginzburg re- 
gion gets absorbed by the Griffiths phase. A complete theory 
of the quantum Griffiths phase will be presented elsewhere. 
We note here that while it is not clear whether the evasive 
Griffiths phase can be observed in real experiment, the quan- 
tum fluctuation effects described by the general equation (fTTT i 
should definitely be accessible at least for 6r > g~ l and possi- 



bly in an even wider range of parameters. 

The physical systems where these fluctuations may be 
experimentally observed include disordered superconducting 
films with an unusual pairing symmetry and possibly high-T c 
cuprates in the vicinity of the termination of the supercon- 
ducting dome in the overdoped regime. Independently of the 
nature of the transition in the overdoped high-T c materials, the 
effective theory still should have the form (JSJ [but the physi- 
cal meaning of the transition parameter r in (|6]l may be differ- 
ent] and therefore all functional dependencies should remain 
the same. We also note here that recent STM experiments, 
(see e.g. Ref.ll7l) suggest that disorder may play an important 
role in high-r c materials and, if so, should contribute to T c - 
suppression [9]. Another superconducting system where the 
field theory of the transition is identical to (|6]l is an s-wave 
superconducting film in which the transition temperature is 
suppressed to zero by magnetic disorder. Finally, we empha- 
size that since a superconductor is a perfect diamagnetic ma- 
terial, fluctuation diamagnetism is physically a correction to 
this perfect diamagnetic behavior and in the vicinity of a su- 
perconducting transition, the magnetic susceptibility is always 
much larger than the Fermi liquid terms. Thus, the experi- 
mental verification of the predicted quantum critical behavior 
of the magnetic susceptibility is certainly feasible. 
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